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Abstract. We study a specific convex maximization problem in n- 
dimensional space. The conjectured solution is proved to be a vertex of the 
polyhedral feasible region, but only a partial proof of local maximality is known. 
Integer sequences with interesting patterns arise in the analysis, owing to the 
number theoretic origin of the problem. 

Dedicated in memory of Jill Spurr Titus, with love 

1. Problem 
For each positive integer n, maximize the convex function 

n 1 

/ rf . 

■ -. J/ J 
1 = 1 

over the polyhedron in n-dimensional real space M n defined by 

(j + l)xj + Xi > (j + l)i + Eij for 1 < j < i < n 
where £jj = 1 if i — j — 1 and e^- = otherwise. Prove that: 

(i) a global maximum (a 1; a 2 , a n ) exists and is unique 

(ii) the components a« of the global maximum satisfy 

cq = 1, a.2 = 2, = 4 

and, when i > 4, 

a { = (j + l)(i - dj) 
for any j with (j + l)aj - ja^-i < i < ( j + 2)a i+1 - (j + l)aj . 

Remark. A solution of this problem will imply the truth of a certain number 
theoretic conjecture due to Levine and O' Sullivan || . 
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2. Partial Solution 
For fixed n, let £ = (xi,x 2 , ...,x n ), 



n 
t=l 

and 

P = {£ ■ {j + l)xj + x l > {j + l)z + eij, l<j<i<n} 

If £ G P, then clearly x^ > 1 for 1 < z < n. As a consequence, P contains no lines 
and / is bounded above on P; therefore, the supremum of / over P is attained at 
one or more vertices of P 0. This proves the existence part of (i). While we do not 
know how to prove the remainder of (i) or (ii), we show here that the conjectured 
global maximum a = (a\, a 2 , a n ) is: 

(a) well-defined 

(b) feasible (that is, a G P) 

(c) a vertex of P 

and, if certain key inequalities hold, 

(d) a local maximum of /(£) subject to £ G P. 

2.1. Proof of (a). The well- definition issue arises because of the conceivable non- 
uniqueness (or even non-existence) of j when determining a, for i > 4. Let 

Cj = (j + 2)a i+1 - (j + l)aj 

for j > 1. Expressed using c\, c%, the definition of a 4 , a 5 ... is 

' 3(i - 2) for 4 = c x < i < c 2 = 10 
4(i - 4) for 10 = c 2 < i < c 3 = 14 
5(i - 6) for 14 = c 3 < z < c 4 = 24 



We prove that both sequences ai, ... and ci, c 2 , ... are strictly increasing. Hypoth- 
esize inductively that c$_i > q_ 2 > ••• > c 2 > c\ = 4, where % > 4 is fixed. Since 
q_i > z", there exists uniquely j < i with Cj_i < z < Cj. If z + 1 < Cj, then 
Cj_i < i + 1 < Cj and hence 

Of+i - at = (j + 1)(* + 1 - a/) - (j + l)(i - aj) = j + 1 > 
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If i + 1 = Cj, then Cj < i + 1 < Cj + \ and hence 

a i+1 -ai = (j + 2)(i + 1 - a j+1 ) - (j + l)(i - aj) 

= [(j + 2)(i + 1) - (j + [(j + 2)a j+1 - (j + l)a 3 ] 

= [(j + 2)(i + l)-(j + l)i]- Cj = [(j + 2)(i + l)-(J + l)i]-(i + l)=j + l>0 

We deduce that (aj+i — a,) — (aj — aj_i) > and thus 

Ci - Ci-i = (i + 2)a i+ i - 2(i + 1)^ + i = 2a i+i + i(a i+1 + aj_i) - 2(i + 1)^ 
> 2a i+1 + 2 i Oj — 2(i + l)aj = 2(a i+1 — aj) > 

This completes the inductive proof, from which well-definition follows immediately. 
As a consequence, we may define 

f 1 if 1 < i < 3 

[ j if A < i < n and Cj_i < i < cj 

without ambiguity. 

2.2. Proof of (b). This is trivial if 1 < % < 3. If i > 4 and j = b { , it follows that 

k-l 

U + !)(«- - 0' + ^ + - = 5^( c ^'+ m ~ *) > 

m=0 

for 1 < k < n — j and 

k 

+ Oj) - (j - A; + l)(i - aj-fe) = - Cj-_ m ) > 

m=l 

for < A; < j — 1. Both series are telescoping and the inequalities are consequences 
of part (a). We deduce that 

a t = (j + l)(i - aj) > (p+l)(i - Op) 

for any 1 < p < n, from which feasibility of a follows immediately. 

2.3. Proof of (c). This is true since a lies at the intersection of the n hyperplanes 

(hi + l)x bi + Xi = (h + l)i + e ibi 1 < i < n 
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2.4. Key Inequalities. Before discussing part (d), we need to state certain key 
inequalities which, although unproven, appear to be true for all n < 10000 via com- 
puter check. 

Definition. Fix integers i and j with 1 < i < j. Let 

h = j-l 
h = b kl 
h = h 2 



k m -l — bk m _ 2 

k m = bk ml 

where m is the smallest integer such that % > k m . Clearly such an integer m exists. 
Then define 



(-l)™JJ(fcp + l) if i = k n 

P =i 

if i > k r , 



Conjecture. Let Co = 2 for convenience, then 

fcn+l / min{ C j_i-l,n} 

x * = ^ + E K E ^1 ^° f( » ; 1 

* j=i+l \ k=Cj- 2 k 

For example, if n = 24, then £* is the solution of the linear system M£ = v , where 
M is the 24 x 24 identity matrix plus some upper triangular entries in the first b n = 5 
rows as indicated: 

( \ 2 2 00000000000000^ 
1 333333 00000000000000 
1 444400000000000 
100000000055555555550 
1 6 
1 

1 
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and v is the 24- vector with i th element 1/af. The inverse, M -1 , of M is given by 
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V 

and the entries x* of = M~ x v are prescribed by the above summation formula. In 
the case n = 24, we compute 

„ _ 123587941503427 
Xl ~ 187646731272000 > 

„ _ 3536905093973 
X2 ~ 27799515744000 > 
„ 44159 
Xs ~ 1016064 > 

, _ 9439261073843 () 
X4 ~ 750586925088000 > 

47 

^ = 4050 >0 

and these positivity results are consistent with the Conjecture. Of course, x* > for 
i > b n immediately. 

2.5. Partial proof of (d). It suffices to solve the following (primal) linear pro- 
gramming problem: 

n 

J. j 



Minimize g(£) = (-£) • V/(a) = ^ 4 

i=l ^ 

subject to £, e Q 

where Q is the polyhedron 

Q = (bi + l)x bi +Xi> (pi + l)i + e ibi , 1 < i < n} 
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Note that Q contains P and possesses a unique vertex, a. Note also that, by the 
Conjecture, the dual linear programming problem has nonempty feasible region 



min{cj — l,n} 

f : Xj + (j + 1) ^2 x i = — for 1 < j < b r , 

a - 



R= < 



l=Cj-! 3 

Xj = for b n < j < n 

a) 

Xj > for all j 



= {M-'v} = {C} 



hence g is bounded below on Q. Therefore a is the global minimum of g(^) subject 
to £ e Q, which implies that a is a local maximum of /(£) subject to £ e P. 

2.6. Partial proof of the Conjecture. The key inequalities are provably true 
when % is sufficiently large relative to n. More precisely, if 

K n <i<n 

then 

Ci — l 



1 3=Ci-i ■> 

To see this, we prove two lemmas. 



Lemma One. d i:i+1 = + for all i > 1 and d {j > if % > b bb and 
i + 1< j < b n + I. 

Proof of Lemma One. The first part is trivial. The second part is proved by 
noting that m > 1 since k\ — j — 1 > i, so either m = 2 (which implies that dij > 0) 
orm = 3 since 

h = b bj _, < b bbn < i 
(which, in turn, implies that dij = 0). QED. 

c,— 1 

Lemma Two. — — ( j + 1) ^ — > for all j > 1. 

Proof of Lemma Two. Direct computation proves the inequality for j = 1,2,3, 4, 10 
and 14. For all other values of j, we will show that 
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that is, 

e 5 = {j + l)c,-_i( Ci _i - 2 a,) + [{j + 1) - ( Cj - Ci _i)] a, 2 > 
which implies the truth of the Lemma. Observe that, if c^-i < j < c k , then 

Cj_i — 2<ij = {k + l)a fc > 

and 



(j + 1) - ( Cj - Cj _x) = | 



j — 2k — 1 > if Cfc_i < j < Cfc 

-2 A; - 3 < if j = Cfc 

These inequalities yield ej > when j ^ Cfc for any k. In the event j = Cfc for some 
k > 4, the argument is only slightly more complicated: 



ej = (k + l) 2 [(c k + l)(2c k -a k )a k -(2k + 3)(c k -a k ) 2 ] 
> (k + l) 2 [(^ - a k ) (2 Cfc - 2 a k )a k - (2 A; + 3) (c k - a k ) 2 ] 
= (k + l) 2 (c k -a k ) 2 [2a k -(2k + 3)}>0 

for all k > 4. QED. 

3. Closing Words 
Techniques for numerical convex maximization abound Q. A vertex enumeration 
scheme has led to verification that a is the global maximum of /(£) subject to £ e P 
for small n only. Keith Briggs has used the general-purpose optimization programs 
AMPL and LANCELOT to confirm the global maximum claim up to n = 24, and 
CFSQP to do likewise up to n = 121. 

The continuous analog of this problem (with summations replaced by integrals) 
is discussed in a companion paper. 



An outcome of Levine and O' Sullivan's work || is that, for any n, there is a global 
maximum a that satisfies a\ = 1, ci2 = 2, 03 = 4 and either CI4 = 6 or CI4 > 28, where £ 
is restricted to integer points in P (that is, to £ e PflZ n ). Their proof unfortunately 
does not extend to the real case. 



Do there exist other functions / and polyhedra P for which the maximizing vertex 
a is " self-generating" as the dimension n increases? A simple characterization of such 
pairs (/, P) may lead to the insight necessary to solve this problem. 
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